Leukemia circulation kinetics revealed through blood exchange method.

Leukemias and their bone marrow microenvironments undergo dynamic changes over the course of disease. However, little is known about the circulation kinetics of leukemia cells, nor the impact of specific factors on the clearance of circulating leukemia cells (CLCs) from the blood. To gain a basic understanding of CLC dynamics over the course of disease progression and therapeutic response, we apply a blood exchange method to mouse models of acute leukemia. We find that CLCs circulate in the blood for 1-2 orders of magnitude longer than solid tumor circulating tumor cells. We further observe that: (i) leukemia presence in the marrow can limit the clearance of CLCs in a model of acute lymphocytic leukemia (ALL), and (ii) CLCs in a model of relapsed acute myeloid leukemia (AML) can clear faster than their untreated counterparts. Our approach can also directly quantify the impact of microenvironmental factors on CLC clearance properties. For example, data from two leukemia models suggest that E-selectin, a vascular adhesion molecule, alters CLC clearance. Our research highlights that clearance rates of CLCs can vary in response to tumor and treatment status and provides a strategy for identifying basic processes and factors that govern the kinetics of circulating cells.

Leukemia circulation kinetics revealed through blood exchange method.

Leukemias and their bone marrow microenvironment are known to undergo dynamic changes over the course of disease. However, relatively little is known about the circulation kinetics of leukemia cells, nor the impact of specific factors on the clearance of circulating leukemia cells (CLCs) from the blood. To gain a basic understanding of leukemia cell dynamics over the course of disease progression and therapeutic response, we apply a blood exchange method to mouse models of acute leukemia. We find that CLCs circulate in the blood for 1-2 orders of magnitude longer than solid tumor circulating tumor cells. We further observe that: i) leukemia presence in the marrow can limit the clearance of CLCs in a model of acute lymphocytic leukemia (ALL), and ii) CLCs in a model of relapsed acute myeloid leukemia (AML) can clear faster than their untreated counterparts. Our approach can also directly quantify the impact of microenvironmental factors on CLC clearance properties. For example, data from two leukemia models suggest that E-selectin, a vascular adhesion molecule, alters CLC clearance. Our research highlights that clearance rates of CLCs can vary in response to tumor and treatment status and provides a strategy for identifying basic processes and factors that govern the kinetics of circulating cells.

Therapy-induced APOBEC3A drives evolution of persistent cancer cells.

Acquired drug resistance to anticancer targeted therapies remains an unsolved clinical problem. Although many drivers of acquired drug resistance have been identified1-4, the underlying molecular mechanisms shaping tumour evolution during treatment are incompletely understood. Genomic profiling of patient tumours has implicated apolipoprotein B messenger RNA editing catalytic polypeptide-like (APOBEC) cytidine deaminases in tumour evolution; however, their role during therapy and the development of acquired drug resistance is undefined. Here we report that lung cancer targeted therapies commonly used in the clinic can induce cytidine deaminase APOBEC3A (A3A), leading to sustained mutagenesis in drug-tolerant cancer cells persisting during therapy. Therapy-induced A3A promotes the formation of double-strand DNA breaks, increasing genomic instability in drug-tolerant persisters. Deletion of A3A reduces APOBEC mutations and structural variations in persister cells and delays the development of drug resistance. APOBEC mutational signatures are enriched in tumours from patients with lung cancer who progressed after extended responses to targeted therapies. This study shows that induction of A3A in response to targeted therapies drives evolution of drug-tolerant persister cells, suggesting that suppression of A3A expression or activity may represent a potential therapeutic strategy in the prevention or delay of acquired resistance to lung cancer targeted therapy.

YAP1 maintains active chromatin state in head and neck squamous cell carcinomas that promotes tumorigenesis through cooperation with BRD4.

Analysis of The Cancer Genome Atlas and other published data of head and neck squamous cell carcinoma (HNSCC) reveals somatic alterations of the Hippo-YAP pathway in approximately 50% of HNSCC. Better strategies to target the YAP1 transcriptional complex are sought. Here, we show that FAT1, an upstream inhibitor of YAP1, is mutated either by missense or by truncating mutation in 29% of HNSCC. Comprehensive proteomic and drug-screening studies across pan-cancer models confirm that FAT1-mutant HNSCC exhibits selective and higher sensitivity to BRD4 inhibition by JQ1. Epigenomic analysis reveals an active chromatin state in FAT1-mutant HNSCC cells that is driven by the YAP/TAZ transcriptional complex through recruitment of BRD4 to deposit active histone marks, thereby maintaining an oncogenic transcriptional state. This study reveals a detailed cooperative mechanism between YAP1 and BRD4 in HNSCC and suggests a specific therapeutic opportunity for the treatment of this subset of head and neck cancer patients.

Extracellular matrix proteins regulate NK cell function in peripheral tissues.

Natural killer (NK) cells reject major histocompatibility complex class I (MHC-I)-deficient bone marrow through direct cytotoxicity but not solid organ transplants devoid of MHC-I. Here, we demonstrate an immediate switch in NK cell function upon exit from the circulation, characterized by a shift from direct cytotoxicity to chemokine/cytokine production. In the skin transplant paradigm, combining an NK cell-specific activating ligand, m157, with missing self MHC-I resulted in complete graft rejection, which was dependent on NK cells as potential helpers and T cells as effectors. Extracellular matrix proteins, collagen I, collagen III, and elastin, blocked NK cell cytotoxicity and promoted their chemokine/cytokine production. NK cell cytotoxicity against MHC-I-deficient melanoma in the skin was markedly increased by blocking tumor collagen deposition. MHC-I down-regulation occurred in solid human cancers but not leukemias, which could be directly targeted by circulating cytotoxic NK cells. Our findings uncover a fundamental mechanism that restricts direct NK cell cytotoxicity in peripheral tissues.

Antibody-Peptide Epitope Conjugates for Personalized Cancer Therapy.

Antibody-peptide epitope conjugates (APEC) are a new class of modified antibody-drug conjugates that redirect T-cell viral immunity against tumor cells. APECs contain a tumor-specific protease cleavage site linked to a patient-specific viral epitope, resulting in presentation of viral epitopes on cancer cells and subsequent recruitment and killing by CD8+ T cells. Here we developed an experimental pipeline to create patient-specific APECs and identified new preclinical therapies for ovarian carcinoma. Using functional assessment of viral peptide antigen responses to common viruses like cytomegalovirus (CMV) in patients with ovarian cancer, a library of 192 APECs with distinct protease cleavage sequences was created using the anti-epithelial cell adhesion molecule (EpCAM) antibody. Each APEC was tested for in vitro cancer cell killing, and top candidates were screened for killing xenograft tumors grown in zebrafish and mice. These preclinical modeling studies identified EpCAM-MMP7-CMV APEC (EpCAM-MC) as a potential new immunotherapy for ovarian carcinoma. Importantly, EpCAM-MC also demonstrated robust T-cell responses in primary ovarian carcinoma patient ascites samples. This work highlights a robust, customizable platform to rapidly develop patient-specific APECs. SIGNIFICANCE: This study develops a high-throughput preclinical platform to identify patient-specific antibody-peptide epitope conjugates that target cancer cells and demonstrates the potential of this immunotherapy approach for treating ovarian carcinoma.

NR4A1 regulates expression of immediate early genes, suppressing replication stress in cancer.

Deregulation of oncogenic signals in cancer triggers replication stress. Immediate early genes (IEGs) are rapidly and transiently expressed following stressful signals, contributing to an integrated response. Here, we find that the orphan nuclear receptor NR4A1 localizes across the gene body and 3' UTR of IEGs, where it inhibits transcriptional elongation by RNA Pol II, generating R-loops and accessible chromatin domains. Acute replication stress causes immediate dissociation of NR4A1 and a burst of transcriptionally poised IEG expression. Ectopic expression of NR4A1 enhances tumorigenesis by breast cancer cells, while its deletion leads to massive chromosomal instability and proliferative failure, driven by deregulated expression of its IEG target, FOS. Approximately half of breast and other primary cancers exhibit accessible chromatin domains at IEG gene bodies, consistent with this stress-regulatory pathway. Cancers that have retained this mechanism in adapting to oncogenic replication stress may be dependent on NR4A1 for their proliferation.

APOBEC Mutagenesis Is Concordant between Tumor and Viral Genomes in HPV-Positive Head and Neck Squamous Cell Carcinoma.

APOBEC is a mutagenic source in human papillomavirus (HPV)-mediated malignancies, including HPV+ oropharyngeal squamous cell carcinoma (HPV + OPSCC), and in HPV genomes. It is unknown why APOBEC mutations predominate in HPV + OPSCC, or if the APOBEC-induced mutations observed in both human cancers and HPV genomes are directly linked. We performed sequencing of host somatic exomes, transcriptomes, and HPV16 genomes from 79 HPV + OPSCC samples, quantifying APOBEC mutational burden and activity in both host and virus. APOBEC was the dominant mutational signature in somatic exomes. In viral genomes, there was a mean of five (range 0-29) mutations per genome. The mean of APOBEC mutations in viral genomes was one (range 0-5). Viral APOBEC mutations, compared to non-APOBEC mutations, were more likely to be low-variant allele fraction mutations, suggesting that APOBEC mutagenesis actively occurrs in viral genomes during infection. HPV16 APOBEC-induced mutation patterns in OPSCC were similar to those previously observed in cervical samples. Paired host and viral analyses revealed that APOBEC-enriched tumor samples had higher viral APOBEC mutation rates (p = 0.028), and APOBEC-associated RNA editing (p = 0.008), supporting the concept that APOBEC mutagenesis in host and viral genomes is directly linked and occurrs during infection. Using paired sequencing of host somatic exomes, transcriptomes, and viral genomes, we demonstrated for the first-time definitive evidence of concordance between tumor and viral APOBEC mutagenesis. This finding provides a missing link connecting APOBEC mutagenesis in host and virus and supports a common mechanism driving APOBEC dysregulation.

RASAL2 Confers Collateral MEK/EGFR Dependency in Chemoresistant Triple-Negative Breast Cancer.

PURPOSE: While chemotherapy remains the standard treatment for triple-negative breast cancer (TNBC), identifying and managing chemoresistant tumors has proven elusive. We sought to discover hallmarks and therapeutically actionable features of refractory TNBC through molecular analysis of primary chemoresistant TNBC specimens. EXPERIMENTAL DESIGN: We performed transcriptional profiling of tumors from a phase II clinical trial of platinum chemotherapy for advanced TNBC (TBCRC-009), revealing a gene expression signature that identified de novo chemorefractory tumors. We then employed pharmacogenomic data mining, proteomic and other molecular studies to define the therapeutic vulnerabilities of these tumors. RESULTS: We reveal the RAS-GTPase-activating protein (RAS-GAP) RASAL2 as an upregulated factor that mediates chemotherapy resistance but also an exquisite collateral sensitivity to combination MAP kinase kinase (MEK1/2) and EGFR inhibitors in TNBC. Mechanistically, RASAL2 GAP activity is required to confer kinase inhibitor sensitivity, as RASAL2-high TNBCs sustain basal RAS activity through suppression of negative feedback regulators SPRY1/2, together with EGFR upregulation. Consequently, RASAL2 expression results in failed feedback compensation upon co-inhibition of MEK1/2 and EGFR that induces synergistic apoptosis in vitro and in vivo. In patients with TNBC, high RASAL2 levels predict clinical chemotherapy response and long-term outcomes, and are associated via direct transcriptional regulation with activated oncogenic Yes-Associated Protein (YAP). Accordingly, chemorefractory patient-derived TNBC models exhibit YAP activation, high RASAL2 expression, and tumor regression in response to MEK/EGFR inhibitor combinations despite well-tolerated intermittent dosing. CONCLUSIONS: These findings identify RASAL2 as a mediator of TNBC chemoresistance that rewires MAPK feedback and cross-talk to confer profound collateral sensitivity to combination MEK1/2 and EGFR inhibitors.

An extended APOBEC3A mutation signature in cancer.

APOBEC mutagenesis, a major driver of cancer evolution, is known for targeting TpC sites in DNA. Recently, we showed that APOBEC3A (A3A) targets DNA hairpin loops. Here, we show that DNA secondary structure is in fact an orthogonal influence on A3A substrate optimality and, surprisingly, can override the TpC sequence preference. VpC (non-TpC) sites in optimal hairpins can outperform TpC sites as mutational hotspots. This expanded understanding of APOBEC mutagenesis illuminates the genomic Twin Paradox, a puzzling pattern of closely spaced mutation hotspots in cancer genomes, in which one is a canonical TpC site but the other is a VpC site, and double mutants are seen only in trans, suggesting a two-hit driver event. Our results clarify this paradox, revealing that both hotspots in these twins are optimal A3A substrates. Our findings reshape the notion of a mutation signature, highlighting the additive roles played by DNA sequence and DNA structure.

Spectrum of Mechanisms of Resistance to Crizotinib and Lorlatinib in ROS1 Fusion-Positive Lung Cancer.

PURPOSE: Current standard initial therapy for advanced, ROS proto-oncogene 1, receptor tyrosine kinase fusion (ROS1)-positive (ROS1+) non-small cell lung cancer (NSCLC) is crizotinib or entrectinib. Lorlatinib, a next-generation anaplastic lymphoma kinase/ROS1 inhibitor, recently demonstrated efficacy in ROS1+ NSCLC, including in crizotinib-pretreated patients. However, mechanisms of lorlatinib resistance in ROS1+ disease remain poorly understood. Here, we assessed mechanisms of resistance to crizotinib and lorlatinib. EXPERIMENTAL DESIGN: Biopsies from patients with ROS1 + NSCLC progressing on crizotinib or lorlatinib were profiled by genetic sequencing. RESULTS: From 55 patients, 47 post-crizotinib and 32 post-lorlatinib biopsies were assessed. Among 42 post-crizotinib and 28 post-lorlatinib biopsies analyzed at distinct timepoints, ROS1 mutations were identified in 38% and 46%, respectively. ROS1 G2032R was the most commonly occurring mutation in approximately one third of cases. Additional ROS1 mutations included D2033N (2.4%) and S1986F (2.4%) post-crizotinib and L2086F (3.6%), G2032R/L2086F (3.6%), G2032R/S1986F/L2086F (3.6%), and S1986F/L2000V (3.6%) post-lorlatinib. Structural modeling predicted ROS1L2086F causes steric interference to lorlatinib, crizotinib, and entrectinib, while it may accommodate cabozantinib. In Ba/F3 models, ROS1L2086F, ROS1G2032R/L2086F, and ROS1S1986F/G2032R/L2086F were refractory to lorlatinib but sensitive to cabozantinib. A patient with disease progression on crizotinib and lorlatinib and ROS1 L2086F received cabozantinib for nearly 11 months with disease control. Among lorlatinib-resistant biopsies, we also identified MET amplification (4%), KRAS G12C (4%), KRAS amplification (4%), NRAS mutation (4%), and MAP2K1 mutation (4%). CONCLUSIONS: ROS1 mutations mediate resistance to crizotinib and lorlatinib in more than one third of cases, underscoring the importance of developing next-generation ROS1 inhibitors with potency against these mutations, including G2032R and L2086F. Continued efforts are needed to elucidate ROS1-independent resistance mechanisms.

HPV+ oropharyngeal squamous cell carcinomas from patients with two tumors display synchrony of viral genomes yet discordant mutational profiles and signatures.

Human papillomavirus (HPV) positive oropharyngeal squamous cell carcinoma (HPV + OPSCC) is increasing in prevalence in the USA, as are cases of patients with multiple HPV + OPSCCs (mHPV + OPSCC). mHPV + OPSCCs present a unique opportunity to examine HPV + OPSCC mutation acquisition and evolution. We performed sequencing of the viral genome, somatic exome and somatic transcriptome from 8 patients each with 2 spatially distinct HPV + OPSCCs, and 37 'traditional' HPV + OPSCCs to first address if paired tumors are caused by the same viral isolate and next, if acquired alterations, and the underlying processes driving mutagenesis, are shared within pairs. All tumor pairs contained viral genomes from the same HPV type 16 sublineage and differed by 0-2 clonal single nucleotide polymorphisms (SNPs), suggesting infection with the same viral isolate. Despite this, there was significant discordance in expression profiles, mutational burden and mutational profiles between tumors in a pair, with only two pairs sharing any overlapping mutations (3/3343 variants). Within tumor pairs there was a striking discrepancy of mutational signatures, exemplified by no paired tumors sharing high APOBEC mutational burden. Here, leveraging mHPV + OPSCCs as a model system to study mutation acquisition in virally mediated tumors, in which the germline, environmental exposures, immune surveillance and tissue/organ type were internally controlled, we demonstrate that despite infection by the same viral isolate, paired mHPV + OPSCCs develop drastically different somatic alterations and even more strikingly, appear to be driven by disparate underlying mutational processes. Thus, despite a common starting point, HPV + OPSCCs evolve through variable mutational processes with resultant stochastic mutational profiles.

Small cell transformation of ROS1 fusion-positive lung cancer resistant to ROS1 inhibition.

Histologic transformation from non-small cell to small cell lung cancer has been reported as a resistance mechanism to targeted therapy in EGFR-mutant and ALK fusion-positive lung cancers. Whether small cell transformation occurs in other oncogene-driven lung cancers remains unknown. Here we analyzed the genomic landscape of two pre-mortem and 11 post-mortem metastatic tumors collected from an advanced, ROS1 fusion-positive lung cancer patient, who had received sequential ROS1 inhibitors. Evidence of small cell transformation was observed in all metastatic sites at autopsy, with inactivation of RB1 and TP53, and loss of ROS1 fusion expression. Whole-exome sequencing revealed minimal mutational and copy number heterogeneity, suggestive of "hard" clonal sweep. Patient-derived models generated from autopsy retained features consistent with small cell lung cancer and demonstrated resistance to ROS1 inhibitors. This case supports small cell transformation as a recurring resistance mechanism, and underscores the importance of elucidating its biology to expand therapeutic opportunities.

Quantification of ongoing APOBEC3A activity in tumor cells by monitoring RNA editing at hotspots.

APOBEC3A is a cytidine deaminase driving mutagenesis, DNA replication stress and DNA damage in cancer cells. While the APOBEC3A-induced vulnerability of cancers offers an opportunity for therapy, APOBEC3A protein and mRNA are difficult to quantify in tumors due to their low abundance. Here, we describe a quantitative and sensitive assay to measure the ongoing activity of APOBEC3A in tumors. Using hotspot RNA mutations identified from APOBEC3A-positive tumors and droplet digital PCR, we develop an assay to quantify the RNA-editing activity of APOBEC3A. This assay is superior to APOBEC3A protein- and mRNA-based assays in predicting the activity of APOBEC3A on DNA. Importantly, we demonstrate that the RNA mutation-based APOBEC3A assay is applicable to clinical samples from cancer patients. Our study presents a strategy to follow the dysregulation of APOBEC3A in tumors, providing opportunities to investigate the role of APOBEC3A in tumor evolution and to target the APOBEC3A-induced vulnerability in therapy.

Aneuploidy and a deregulated DNA damage response suggest haploinsufficiency in breast tissues of BRCA2 mutation carriers.

Women harboring heterozygous germline mutations of BRCA2 have a 50 to 80% risk of developing breast cancer, yet the pathogenesis of these cancers is poorly understood. To reveal early steps in BRCA2-associated carcinogenesis, we analyzed sorted cell populations from freshly-isolated, non-cancerous breast tissues of BRCA2 mutation carriers and matched controls. Single-cell whole-genome sequencing demonstrates that >25% of BRCA2 carrier (BRCA2mut/+ ) luminal progenitor (LP) cells exhibit sub-chromosomal copy number variations, which are rarely observed in non-carriers. Correspondingly, primary BRCA2mut/+ breast epithelia exhibit DNA damage together with attenuated replication checkpoint and apoptotic responses, and an age-associated expansion of the LP compartment. We provide evidence that these phenotypes do not require loss of the wild-type BRCA2 allele. Collectively, our findings suggest that BRCA2 haploinsufficiency and associated DNA damage precede histologic abnormalities in vivo. Using these hallmarks of cancer predisposition will yield unanticipated opportunities for improved risk assessment and prevention strategies in high-risk patients.

MET Alterations Are a Recurring and Actionable Resistance Mechanism in ALK-Positive Lung Cancer.

PURPOSE: Most ALK-positive lung cancers will develop ALK-independent resistance after treatment with next-generation ALK inhibitors. MET amplification has been described in patients progressing on ALK inhibitors, but frequency of this event has not been comprehensively assessed. EXPERIMENTAL DESIGN: We performed FISH and/or next-generation sequencing on 207 posttreatment tissue (n = 101) or plasma (n = 106) specimens from patients with ALK-positive lung cancer to detect MET genetic alterations. We evaluated ALK inhibitor sensitivity in cell lines with MET alterations and assessed antitumor activity of ALK/MET blockade in ALK-positive cell lines and 2 patients with MET-driven resistance. RESULTS: MET amplification was detected in 15% of tumor biopsies from patients relapsing on next-generation ALK inhibitors, including 12% and 22% of biopsies from patients progressing on second-generation inhibitors or lorlatinib, respectively. Patients treated with a second-generation ALK inhibitor in the first-line setting were more likely to develop MET amplification than those who had received next-generation ALK inhibitors after crizotinib (P = 0.019). Two tumor specimens harbored an identical ST7-MET rearrangement, one of which had concurrent MET amplification. Expressing ST7-MET in the sensitive H3122 ALK-positive cell line induced resistance to ALK inhibitors that was reversed with dual ALK/MET inhibition. MET inhibition resensitized a patient-derived cell line harboring both ST7-MET and MET amplification to ALK inhibitors. Two patients with ALK-positive lung cancer and acquired MET alterations achieved rapid responses to ALK/MET combination therapy. CONCLUSIONS: Treatment with next-generation ALK inhibitors, particularly in the first-line setting, may lead to MET-driven resistance. Patients with acquired MET alterations may derive clinical benefit from therapies that target both ALK and MET.

A post-transcriptional program of chemoresistance by AU-rich elements and TTP in quiescent leukemic cells.

BACKGROUND: Quiescence (G0) is a transient, cell cycle-arrested state. By entering G0, cancer cells survive unfavorable conditions such as chemotherapy and cause relapse. While G0 cells have been studied at the transcriptome level, how post-transcriptional regulation contributes to their chemoresistance remains unknown. RESULTS: We induce chemoresistant and G0 leukemic cells by serum starvation or chemotherapy treatment. To study post-transcriptional regulation in G0 leukemic cells, we systematically analyzed their transcriptome, translatome, and proteome. We find that our resistant G0 cells recapitulate gene expression profiles of in vivo chemoresistant leukemic and G0 models. In G0 cells, canonical translation initiation is inhibited; yet we find that inflammatory genes are highly translated, indicating alternative post-transcriptional regulation. Importantly, AU-rich elements (AREs) are significantly enriched in the upregulated G0 translatome and transcriptome. Mechanistically, we find the stress-responsive p38 MAPK-MK2 signaling pathway stabilizes ARE mRNAs by phosphorylation and inactivation of mRNA decay factor, Tristetraprolin (TTP) in G0. This permits expression of ARE mRNAs that promote chemoresistance. Conversely, inhibition of TTP phosphorylation by p38 MAPK inhibitors and non-phosphorylatable TTP mutant decreases ARE-bearing TNFα and DUSP1 mRNAs and sensitizes leukemic cells to chemotherapy. Furthermore, co-inhibiting p38 MAPK and TNFα prior to or along with chemotherapy substantially reduces chemoresistance in primary leukemic cells ex vivo and in vivo. CONCLUSIONS: These studies uncover post-transcriptional regulation underlying chemoresistance in leukemia. Our data reveal the p38 MAPK-MK2-TTP axis as a key regulator of expression of ARE-bearing mRNAs that promote chemoresistance. By disrupting this pathway, we develop an effective combination therapy against chemosurvival.

Passenger hotspot mutations in cancer driven by APOBEC3A and mesoscale genomic features.

Cancer drivers require statistical modeling to distinguish them from passenger events, which accumulate during tumorigenesis but provide no fitness advantage to cancer cells. The discovery of driver genes and mutations relies on the assumption that exact positional recurrence is unlikely by chance; thus, the precise sharing of mutations across patients identifies drivers. Examining the mutation landscape in cancer genomes, we found that many recurrent cancer mutations previously designated as drivers are likely passengers. Our integrated bioinformatic and biochemical analyses revealed that these passenger hotspot mutations arise from the preference of APOBEC3A, a cytidine deaminase, for DNA stem-loops. Conversely, recurrent APOBEC-signature mutations not in stem-loops are enriched in well-characterized driver genes and may predict new drivers. This demonstrates that mesoscale genomic features need to be integrated into computational models aimed at identifying mutations linked to diseases.

Stromal Microenvironment Shapes the Intratumoral Architecture of Pancreatic Cancer.

Single-cell technologies have described heterogeneity across tissues, but the spatial distribution and forces that drive single-cell phenotypes have not been well defined. Combining single-cell RNA and protein analytics in studying the role of stromal cancer-associated fibroblasts (CAFs) in modulating heterogeneity in pancreatic cancer (pancreatic ductal adenocarcinoma [PDAC]) model systems, we have identified significant single-cell population shifts toward invasive epithelial-to-mesenchymal transition (EMT) and proliferative (PRO) phenotypes linked with mitogen-activated protein kinase (MAPK) and signal transducer and activator of transcription 3 (STAT3) signaling. Using high-content digital imaging of RNA in situ hybridization in 195 PDAC tumors, we quantified these EMT and PRO subpopulations in 319,626 individual cancer cells that can be classified within the context of distinct tumor gland "units." Tumor gland typing provided an additional layer of intratumoral heterogeneity that was associated with differences in stromal abundance and clinical outcomes. This demonstrates the impact of the stroma in shaping tumor architecture by altering inherent patterns of tumor glands in human PDAC.